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I.  lOTRODUCTION 

The  recent  conflict  In  Southeast  Asia  dramatically  revealed  the  vult.era- 
blllty  of  sophisticated  and  expensive  aircraft  to  small  arms  ground  fire  and 
to  surface-to-air  missiles.  Since  the  fuel  tanks  of  tactical  aircraft  have  the 
largest  presented  area  of  all  of  the  vulnerable  components,  engineering  estimates 
of  fuel  tank  response  to  penetrating  ballistic  projectiles  are  required  in  order 
to  design  more  survivable  tanks.  Projectiles  that  penetrate  tanks  containing 
fuel  caiise  more  severe  damage  to  the  tank  than  those  that  pass  through  empty 
tanks.  The  reason  for  this  Is  that  the  passage  of  the  projectile  through  the 
fuel  causes  an  Intense  pressure  pulse  to  propagate  In  the  fuel  and  strike  the 
walls  of  the  tank.  This  large  Internal  fluid  pressure  on  the  walls  can  cause 
severe  petallng  of  the  walls,  usually  at  the  entrance,  and  exit  points  of  the 
projectiles.  The  development  of  the  pressure  pulse  in  the  fuel  by  the  ballistic 
penetrator  Is  known  as  the  hydraulic  ram  effect,  and  the  fluid  pressure  is 
referred  to  as  the  hydraulic  ram  pressure. 

The  hydraulic  ram  effect  can  be  divided  into  three  phases;  the  early  shock 
phase,  the  later  drag  phase,  and  the  cavity  phase.  The  shock  phase  Is 
Initiated  when  the  pr&jectlle  penetrates  the  wall  and  Impacts  the  fluid.  As 
energy  Is  transferred  to  the  fluid,  a strong  hemispherical  shock  wave  centered 
at  the  point  of  impact  Is  formed.  This  creates  an  impulsive  load  on  the  Inside 
of  the  entry  wall  in  the  vicinity  of  the  entry  hole  which  may  cause  the  entry 
wall  to  crack  and  petal.  As  the  projectile  travels  through  the  fluid.  Its 
energy  is  transformed  into  kinetic  energy  of  fluid  motion  as  the  projectile  is 
slowed  by  viscous  drag.  A pressure  field  Is  generated  as  fluid  is  displaced 
from  the  projectile  path.  In  contrast  to  the  pressures  developed  In  the  shock 


pha9e,  the  fluicj  is  accelerated  gradually  rather  than  impulalvelyi  so  that  the 
peak  pressure  Is  much  lower;  however,  the  duration  of  the  pressure  pixlse  is 
considerably  longer.  A cavity  develops  behind  the  projectile  as  it  passes 
ti^rough  the  fluid  which  is  filled  with  liquid  vapor  evaporated  from  the  cavity 
surface  and  air  which  can  enter  the  cavity  through  the  entry  hole.  As  the  fluid 
seeks  to  regain  its  undiscnrbed  condition,  the  cavity  will  oscillate.  The  con-* 
comitant  pressures  will  pump  fluid  from  any  holes  in  the  tank  and  they  may  he 
sufficient  to  damage  fuel  cell  components.  This  cavity  oscillation  is  called 
the  "cavity  phase." 

The  structural  response  of  the  fuel  tank  walls  to  the  hydraulic  ram 
pressure  is  a complicated  process.  The  pressure  in  the  fuel  caused  by  the 
penetrating  projectile  acts  on  the  tank  walls,  causing  them  to  displace.  This 
displacement  in  turn  affects  the  pressure  in  the  fuel,  thus  leading  to  a complex 
interaction  between  the  fuel  and  the  tank  walls.*  This  interaction  phenomenon 
is  referred  to  as  fluid-structure  interaction. 

Considerable  research  effort  has  been  expended  to  develop  an  understanding 
of  the  hydraulic  ram  effect  and  the  associated  fluid-structure  interaction. 
References  1 and  2 contain  extensive  bibliographies  of  most  of  the  work  on 
hydraulic  ram.  Of  particular  interest  here  is  the  Hydraulic  Ram  Project  that 
has  been  conducted  at  the  Naval  Weapons  Center  (NWC) , China  Lake,  California, 

The  first  phase  of  that  project  involved  the  characterization  of  the  pressure 
waves,  generated  by  the  projectile  during  fluid  penetration  through  a rectangular 
tank.  An  analytical  mo iel  for  predicting  the  drag  phase  pressure  field  was 
presented  in  Ref.  3,  During  FY  1973,  an  extens:f  'e  series  of  ballistic  tests 
were  performed  at  NWC  to  obtain  detailed  fluid  pressure  measurements  at 

*Any  cracking  and  petaling  of  the  walls  will  also  change  the  pressure  in  the 
fluid,  and  hence  the  subsequent  loading  on  the  walls. 
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were  performed  at  NWC  to  obtain  detailed  fluid  pressure  measurements  at 

several  locations  for  a variety  of  projectiles  under  a wide  range  of  impact 

* 

conditions.  The  results  of  these  testa  is  reported  in  Ref.  4.  A digital 
computer  code  for  predicting  the  drag  phase  fluid  pressure  in  a rectangular 
tank  due  to  ballistic  penetrators,  based  upon  the  theory  of  Ref.  3 and  the 
empirical  data  of  Ref.  4,  was  subsequently  developed  by  Lundstrom  and  Fxmg. 
Informa^tion  on  the  code  and  its  use  is  given  in  Ref.  5.  This  code  will  be 
referred  to  hereafter  as  the  NWC  hydraulic  ram  pressure  code.  One  very  signifi- 
cant feature  of  this  code  is  the  fact  that  the  wall  surfaces  of  the  rectangular 
fluid  voluioe  are  treated  either  as  free  surfaces  or  as  fixed  surfaces  when 
considering  wave  reflections  from  the  walls. 

The  second  phase  of  the  NWC  Hydraulic  Ram  Project  involved  the  characteri- 
zation of  the  fluid-structure  interaction  during  the  loading  of  the  fuel  tank 
walls  by  the  hydraulic  ram  pressure.  In  support  of  this  goal,  an  extensive 
series  of  ballistic  penetration  tests  were  conducted  on  fluid-filled,  rectangular 
tanks  in  March,  1974,  to  obtain  detailed  measurements  of  the  strains  at  several 
locations  on  the  exit  wall  of  the  tank.  A detailed  description  of  the  test 
set-up  and  the  complete  set  of  results  are  given  in  Ref.  6. 

An  analytical  and  experimental  hydraulic  ram  program  has  been  conducted  at 
the  Naval  Postgraduate  School  (NPS)  in  conjunction  with  the  NWC  program  for 
approximately  three  years.  A ballistic  range  has  been  built  consisting  of  a 
22  caliber  rifle  and  a cubic  fluid-containing  tank.  Fluid  pressures  and  entry 
wall  strains  have  been  measured  for  various  projectile  sizes  and  energy  levels. 
The  results  of  these  experiments  are  given  in  Refs.  7-17.  The  analytical  phase 
of  the  NPS  program  has  been  devoted  to  the  development  of  methods  of  analysis 
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that  will  accurately  predict  the  structural  response  of  the  ta^.dc  walls  to  the 
hydraulic  ran  loading.  The  results  of  this  effort  are  given  in  Refs.  10  and 
17-21. 

As  a consequence  of  the  complexity  of  the  fluid-structure  Interaction  while 
the  projectile  Is  penetrating  the  tank,  and  after  It  exits  the  tank,  an  approxi- 
mate theory  was  proposed  In  Ref.  10  that  allows  the  solution  for  the  fluid 
pressures  to  be  obtained  separately  from  the  solution  for  the  wall  response.* 
This  theory  is  known  as  the  "piston  theory."  It.  Is  an  exact  theory  for  the 
one-dlmenslonal  acoustic  fluid -structure  Interaction  problem.  In  the  piston 
theory  the  response  of  the  wall  Is  computed  using  the  conventional  structural 
response  equations,  with  the  normal  pressure  on  the  wall  q given  by  2p^-pcw 
where  p^  Is  the  Incident  fluid  pressure  at  the  wall,  p is  the  fluid  density, 
c is  the  acoustic  velocity  in  the  fluid,  and  w Is  the  wall  velocity**.  The 
values  for  p^  during  the  drag  phase  can  be  obtained  from  the  NWC  hydraulic 
ram  pressure  code.  For  the  shock  phase  pressure,  the  theory  developed  by 
YurkoA'lch,  (Ref.  22,)  can  be  tised  to  determine  the  incident  pressure  field. 

Two  structural  response  digital  computer  codes,  BR-1  and  SATANS,  have  been 
modified  by  the  author  to  Include  the  piston  theory.  The  modification  consisted 
of  incorporating  the  -pew  loading  component  into  the  equations  of  motion 
of  the  structure.  The  BR-1  code.  Ref.  23,  is  a finite  element  code  for  pre- 
dicting the  inelastic,  large  deflection,  transient  response  of  combat  aircraft 
skin-rib-stringer  structures  when  subjected  to  internal  alrblast  loading.  The 
finite  elements  are  flat  rectangular  plates  and  beam  stiffeners.  The  modified 

*An  exact  method  for  determining  the  response  of  the  fluid-structural  system 
was  presented  in  Ref.  3.  However,  this  method  requires  simultaneous  solution 
for  the  fluid  pressures  and  the  wall  deflections.  This  approach  does  not 
appear  to  be  practical  at  this  time. 

**A  dot  above  a variable  denotes  a derivative  with  respect  to  time. 
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code  Is  called  BR-IHR.  A description  of  the  modification  to  BR-1  and  a User's 
Manual  for  BR-IHR  are  presented  in  Ref.  24.  The  code  SATANS  is  a finite 
difference  - Fourier  series  code  for  the  geometrically  nonlinear  static  and 
dynamic  analysis  of  arbitrarily  loaded  shells  of  revolution.  A detailed 
description  of  SATANS  is  given  in  Ref.  25.  Instructions  for  its  use  are  given 
in  Ref.  26.  For  this  studj'  the  code  is  used  to  analyze  flat  circular  plates. 

This  report  presents  the  results  of  our  efforts  to  determine  the  accuracy 
of  the  structural  response  codes  with  piston  theory  when  predicting  the  entry 
and  exit  wall  strains  due  to  the  fluid  pressures  predicted  by  the  NWC  hydraulic 
ram  pressure  code  and  the  Yurkovich  shock  pressure  theory,  A modification  to 
the  piston  theory  is  also  considered  that  accounts  for  cavitation  at  the 
Interface  due  to  tension  in  the  fluid.  Section  II  contains  a description  of 
the  NWC  exit  wall  and  the  NFS  entry  wall  test  programs.  Section  III  presents 
various  aspects  of  the  fluid-structure  interaction,  including  cavitation  at  the 
interface,  and  Section  IV  describes  the  features  of  BR-lHR  and  SATANS  and  the 
proposed  structural  models  of  the  tank  walls.  Section  V describes  the  NWC 
hydraulic  ram  pressure  code  and  presents  the  predicted  pressure  on  the  entry 
and  exit  walls.  The  entry  and  exit  wall  strains  computed  using  SATANS  with 
piston  theory  and  piston  theory  with  cavitation  are  presented  and  compared 
to  the  measured  strains  in  Section  VI,  The  conclusions  with  regard  to  the 
accuracy  of  the  analyses  are  given  in  Section  VII. 
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II.  EXPERIMENTAL  PROGRAMS 


A.  Naval  Weapons  Center  Exit  Wall  Test  Program 

Mr.  Wallace  Fung  (Code  5114)  has  conducted  a series  of  hydraulic  ram  effects 
tests  on  simulated  aircraft  fuel  tanks  at  the  Naval  Weapons  Center  (NWC),  China 
Lake,  California,  In  which  12.7mm  API  ballistic  projectiles  were  fired  at  a 
fluid  filled,  rectangular  tank.  The  tank  is  60  inches  wide  by  60  Inches  high, 
with  22  inches  between  the  entry  and  exit  walls.  The  top  of  the  tank  Is  open. 

A schema  of  the  test  set-up  Is  shown  In  Fig.  1.  The  central  portion  of  the  entry 
wall  Is  a stretched  rubber  membrane.  The  central  portion  of  the  exit  wall  Is  a 
20  inch  by  20  Inch  by  0.125  Inch  2024-T3  aluminum  test  plate  clamped  to  a 0,25 
inch  steel  wall.  The  edge  clamping  was  obtained  by  compression  (friction) 
between  two  rubber  gaskets  around  the  outer  perimeter  of  the  plate.  Additional 
edge  support  was  provided  by  a pin  through  the  plate  at  two  diametrically 
opposed  corners. 

Fourteen  12,7mm  API  projectiles  were  fired  at  the  tank  with  velocities 
between  1300  fps  and  2900  fps.  The  exit  point  of  each  shot  is  shown  in  Fig,  2, 
The  strains  on  the  wet  and  dry  surfaces  of  the  exit  wall  were  measured  at  the 
four  locations  shown  in  Fig.  2 for  each  shot.  The  strains  at  each  location 
were  measured  in  the  two  directions  Indicated  by  the  short  lines  in  Fig.  2. 
Pressures  were  measured  at  five  locations  along  the  trajectory  and  high-speed 
motion  pictures  were  taken  looking  at  the  dry  side  of  the  exit  wall.  Table  1 
lists  the  projectile  velocity  and  the  extent  of  damage  for  each  shot.  In 
general,  low  velocity  projectiles  caused  some  bulging  of  the  exit  wall 
accompanied  by  short  cracks  approximately  three  inches  in  length.  The  high 
velocity  projectiles  caused  severe  tearing  of  the  wall  with  several  cracks 
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Exit  Points  and  the  Strain  Gages  on  the 


Shot 

Velocity  (fps) 

1 

2,200 

2 

2,217 

3 

2,254 

4 

2,198 

5 

2,113 

6 

2,763 

7 

2,673 

8 

2,643 

9 

2,718 

10 

1,347 

11 

1,386 

12 

1,366 

13 

1,396 

14 

1,389 

Velocity  and 


clt  Wall  Damage 


Damage 

Cigar-shaped  hole  3 in.  long 
Small  circular  hole  2 a.n.  in  dla. 
Cigar-shaped  hole  with  one  19  in.  crack 
Clgar-s>.iaped  hole  2 in.  long 

II  II  II  It  II  II 

Circular  hole  3 in.  in  dla  with  two 
10  in.  cracks 

5 In.  hole  with  several  5 in.  cracks 
and  petallng 

5 In.  hole  with  several  5 in.  cracks 
and  petallng 

Bulged  plate  with  3 in.  cigar-shaped 
hole 

Cigar-shaped  hole  3 in.  long 

II  II  II  II  II  II 

II  II  II  II  II  It 

Not  available 

Very  small  hole  with  two  2 in.  cracks 
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running  from  the  exit  hole  to  the  edges  of  the  plate.  Figure  3,  taken  from 
Ref.  6,  shows  the  test  panel  from  shot  14.  A more  detailed  description  of  the 
experimental  set-up,  a presentation  of  all  of  the  exit  wall  strain  data,  and 
photographs  of  each  of  the  damaged  plates  are  given  In  Ref.  6.  Selected  re- 
sults from  the  NWC  tests  will  be  presented  In  Section  VI  where  the  comparison 
between  the  experimental  rssults  and  the  analytical  results  Is  made. 


B.  Naval  Postgraduate  School  Entry  Wall  Test  Program 

A series  of  experiments  have  been  conducted  at  NPS  using  a simple  tank 
configuration  to  determJ.ne  hydraulic  ram  pressure  loadings  near  the  entry  wall 
for  several  .222  callbc:r  energy  levels.  The  test  fuel  cell  is  cubical,  with 
inside  dimensions  of  17  inches  on  a side,  and  W6S  constructed  by  welding  one 
quarter  inch  by  three  inch  aluminum  angle  sections  to  form  a frame.  Side  and 
bottom  walls  are  one  inch  thick  plexiglass  and  the  entry  and  exit  walls  \;ere 
belted  to  the  frame  so  that  various  thicknesses  and  materials  could  be  tested. 

The  top  is  open  and  bullet  entry  into  the  tank  was  accomplished  by  providing  a 
one  inch  diameter  hole  in  the  center  of  the  entry  wall.  A .222  caliber  rifle 
was  used  to  fire  copper-jacketed  lead  projectiles  whose  mass  and  impact 
velocity  were  adjusted  to  obtain  the  various  impact  energy  levels  tested. 

Spark  shadow(;raphs  were  taken  to  determine  shock  phase  wave  shape,  shock  posi- 
tion and  proj(»ctile  location  as  a function  of  time,  Ref.  13,  Kistler  603H 
quartz  presLure  transducers  were  mounted  in  the  tank  to  measure  internal 
pressure  at  a given  spatial  position  verses  time,  pressures  were  measured 
along  a line  inclined  50®  from  the  vc  tical  at  radii  varying  from  two  to  eight 
inches  from  the  impact  point.  The  internal  pressures  were  measured  using  both 
0.5  inch  and  0.05  inch  thick  steel  plates  at  the  entry  wall  In  order  to  stimulate 
the  case  of  a rigid  entry  wall  and  a flexible  entry  wall  respectively.  The 
pressure  measurements  are  compared  with  the  pressure  predicted  by  the  IWC 
hydraulic  ram  pressure  code  in  Ref.  17.  The  projectile  position  data  given  in 
Ref.  13  was  used  to  determine  the  drag  coefficients  used  In  the  NWC  code. 
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A series  of  e:.perimen£s  to  determine  the  entry  wall  structural  response 
to  the  shock  and  drag  phase  loadings  for  two  energy  levels  (7,493  In, -lb  and 
12,323  In. -lb)  have  also  been  conducted.  The  entry  wall  was  a circular  plate 
with  a fifteen  Inch  diameter  and  a one  Inch  diameter  hole  at  the  center)^.  The 
test  plate  was  bolted  to  a one  Inch  thick  aluminum  plate  with  a fifteen  Inch 
circular  cutout.  7075-T6  Aluminum  plates  0.05,  0.09,  and  0.16  Inches  thick 
were  tested.  Each  entry  wall  had  four  EA-13  type  175  ohm  strain  gages  mounted 
2 inches  from  the  impact  point.  One  pair  of  gages  measured  radial  strain  with 
one  gage  on  the  Insjde  and  the  other  on  the  outside  of  the  plate,  A secono 
pair  of  gages  was  oriented  circumferentially.  Figure  4 shows  the  test  tank 
set-up.  A simple  constant  current  strain  gage  circuit  was  used  to  measure 
the  dynamic  response  of  the  entry  wall.  The  recording  system  consisted  of  a 
dual-beam  oscilloscope  with  two  1A7  pre-amp  plug-in  units  with  differential 
Input  capability  and  a Polaroid  scope  camera.  The  details  of  the  test  and 
thr  oscilloscope  pictures  of  the  strain  traces  are  given  in  Refs.  15  and  17. 
The  shots  at  the  high  energy  level  caused  considerable  tank  leakage,  strain 
gage  bond  failure  and  permanent  deformation  of  the  0,05  In.  plate.  The  low 
energy  shots  did  not  cauf  i any  apparent  damage.  A preliminary  comparison  of 
the  measured  strains  with  the  strains  predicted  by  SATANS  is  made  in  Refs. 

15  and  x7.  A more  detailed  comparison,  including  the  effects  of  cavitation, 
is  made  in  Section  VI  of  this  report. 


i 
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*A  circular  plate  was  tested  so  that  the  results  from  SATANS  are  directly 
applicable. 
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Ill  FLUID-STRUCTURE  INTERACTION 


A.  Governing  Equations 

Consider  the  fluid-structure  interface  at  a typical  portion  of  the 
entrance  and  exit  walls  of  a fuel  tank  shown  in  Fig,  5.  The.  linear 
differential  equation  representing  unsteady,  low  velocity,  irrotational, 
frictionless,  compressible  fluid  flow  can  be  given  in  terms  of  a velocity 
potential  ^ in  the  form 


+ <j)  +4*  " “TO  •I’ 

^xx  ^yy  ^zz  ^2  ^ 


(3) 


where  c Is  the  sonic  velocity  of  the  fluid  (assumed  to  he  constant) , subscripts 


denote  partial  derivatives,  anJ  (j)^,  and  <|>^  are  the  velocities  in  the  x, 


y and  z directions  respectively.  The  pressure  in  the  fluid  p can  he  obtained 
using  Bernoulli's  equation 


p _ _ p|  - £ 


(2) 


3 


where  p is  the  fluid  density,  assumed  constant, 

2 2 2 2 

+ (j)  ^ (3) 

X y z 

and  p is  the  ambient,  pressure,  p is  positive  in  compression.  For  low 

A 

2 

velocity  flow  the  V term  in  Eq.  (2)  is  usually  neglected,  giving 

p _ p « _ (4) 

O 

At  the  wall,’  the  Interface  conditions  are 


p * ± q 


(5) 
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FIG.  5 Typical  portion  of  Entrance  and  Exit  Walls 
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where  q Is  the  normal  pressure  acting  on  the  wall,  and  w Is  the  wall 
deflection  normal  to  the  plane  of  the  wall.  The  minus  sign  on  q applies  at 
the  entry  wall  and  the  plus  sign  applies  at  the  exit  wall.  These  Interface 
''ondltlons  are  based  upon  the  assumption  that  both  compression  and  tension 
pressures  can  exist  at  the  fluid-structure  ln^..rface.  In  reality,  cavitation 
may  occur  and  the  wall  velocity  and  fluid  velocity  may  not, be  the  same. 

A procedure  for  determining  the  pressure  on  the  tank  walls  due  to  transient 
wall  motion,  based  upon  the  above  theory,  was  presented  In  Ref.  3.  This  proce- 
dure used  Klrchoff's  relation  for  the  p .essure  at  any  point  within  the  fluid 
due  to  motion  of  the  surface  of  the  fluid  volume.  The  method  requires  simul- 
taneous solution  for  the  fluid  pressures  and  wall  deflections.  It  was 
emphasized  in  Ref.  3 that  the  proposed  method,  although  simple  In  concept,  Is 
extremely  difficult  to  carry  out  in  practice,  due  mainly  to  computer  storage 
space  and  running  time. 
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B,  One  - Dinensional  Interaction  (Piston  Theory) 

Due  to  the  considerable  difficulties  In  treating  the  fully  coupled 
fluid-structure  Interaction  problem,  the  use  of  simplified  procedures  Is 
appealing.  One  such  procedure  Is  the  piston  theory.  This  theory  has  been 
in  use  since  the  early  1940* s when  It  was  applied  to  the  study  of  the  effect 
of  underwater  explosions  on  ship  plates.  Ref.  27.  It  provides  the  correct  inter- 
action solution  to  the  one-dimensional  propagation  of  plane  stress  waves  in  an 
acoustic  medium  due  to  a moving  plane  botmdary.  Several  recent  studies  have 
been  made  to  determine  Its  accuracy  when  applied  to  two-dimensional  fluid- 
structure  interaction  problems,  such  as  the  propagation  of  plane  waves  over 
circular  cylindrical  shells  and  the  vibration  of  spherical  shells.  Refs,  28 
and  29. 

A derivation  of  the  piston  theory  is  given  In  Kef.  10.  The  essential 
feature  of  piston  theory  is  that  the  normal  pressure  at  the  wall  Is  given  by 

<1  ■ Pq  pc(v^  - w)  (6) 

where  p^  and  v^  are  the  Incident  pressure  and  normal  velocity  of  the 
fluid  at  the  wall  respectively,  and  w Is  the  wall  velocity*.  The  pressure 
p^  and  velocity  v^^  are  the  pressure  and  velocity  that  would  exist  In  the 
fluid  if  the  wall  was  not  there,  l.e,  p^  and  v^  do  not  contain  any  reflec- 
tion effects  due  to  the  presence  of  the  wall.  However,  effects  on  p^ 
and  v^  due  to  reflections  from  other  surfaces  should  be  Included. 

*The  derivation  in  Ref.  10  used  the  relationship  p ■=  - p<p.  If  Eq.  (2)  is 

. 2.  . 2 . 

used,  Eq.  (6)  becomes  q “ p^  + pc(v^  - w)  - pw  . However,  since  w « c, 
the  velocity  squared  term  can  be  neglected. 
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The  pressure  given  by  Eq.  (6)  can  be  interpreted  as  being  composed  of 
two  parts,  p^  , the  incident  pressure  in  thj^  fluid,  and  pc(v^  - w) , the 
reflected  pressure.  Note  that  the  reflected  pressure  is  directly  proportional 
to  (v^  - w) , the  mismatch  between  the  fluid  velocity  and  the  wall  velocity. 

The  pev^  term  is  the  reflected  pressure  in  a plane  wave  at  a rigid,  non-moving 


wall  (w  ■=  0)  and  can  be  shown  to  be  equal  to  p^  , the  incident  pressure.  Thus, 


Eq.  (6)  can  also  be  given  in  the  form 


q - 2p^  - pew 


(7) 


If  the  plane  wave  front  is  moving  normal  to  a rigid,  non-moving  wall,  as 
is  the  case  for  a hemispherlcally  expanding  shock  wave  at  the  entry  wall,  then 
there  is  no  mismatch  between  the  normal  fluid  velocity  and  the  wall  velocity, 


and  the  incident  pressure  p^  should  be  used  in  place  of  the  2p^  term  in 


Eq.  (7). 
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C,  Sprerlcal  Wave  Propagation  and  Interaction  with  a Planar  Wall 

The  NWC  hydraulic  rejn  pressure  code  docs  not  assume  plane  wave  stress 
propagation  in  the  fluid,  but  Instead  treats  the  problem  of  the  penetrating 
projectile  as  a seq^ience  of  spherical  stress  waves  In  an  acoustic  medium  whose 
tenters  of  disturbance  lie  along  the  bullet  path.  Thus,  spherical  waves 
strike  the  planar  tank  wall,  not  planar  waves.  The  pressure  at  a rigid  wall, 
corresponding  to  the  2p^  pressure  In  the  one-dimensional  interaction,  can  he 
determined  using  the  method  of  Images.  If  the  Image  source  Is  of  the  same 
sign  as  the  actual  source,  then  the  normal  velocity  at  the  wall  is  zero,  and 
the  pressure  at  the  wall  Is  the  net  result  of  the  two  sources.  Thus.,  If  Fq.« 
(A)  is  used  for  the  pressure,  the  pressure  doubles  at  the  rigid  wall.  Just  as 
in  the  one-dlmenslonal  planar  wave  case,  and  2p^  is  the  pressuie  at  the 
rigid  wall.  If  the  image  source  Is  of  the  opposite  sign,  the  pressure  Is  ze.:'' 


at  the  wall  and  the  wall  is  like  a free  surface. 


D.  Cavitation 

The  assumption  has  been  made  in  the  discussion  above  that  there  is 
continuity  between  the  fluid  and  the  vjall  at  the  interface.  This  assumption 
becomes  questionable  when  tension  exists  at  the  interface,  l.e.,  when 

p “ 2p^  - pew  < 0 (8) 

Experiments  have  shown  that  cavitation  occurs  in  the  water  in  front  of 
thin,  circular,  air-backed  plates  when  they  are  subjected  to  underwater  ex- 
plosive-type shock  waves.  Ref.  27.  Two  values  obtained  for  the  maximum 
tension  in  seawater  are  70  and  40  psi.  Ref.  30.  Cole  has  suggested  that  it 
usually  is  a good  approximation  to  assume  that  cavitation  begins  at  the 
surface  of  a free  plate,  neglecting  diffraction  effects  from  the  plate  edge, 
when  the  interface  pressure  becomes  zero.  Ref.  30. 


IV  STRUCTURAL  RESPONSE  COIiPUTER  CODES  AND  MODELS 


A.  BR-IHR 

The  Northrup  Corporation,  under  Air  Force  funding,  has  developed  a finite 
element  digital  computer  code,  called  BK-1,  for  predicting  the  inelastic,  large 
deflection,  transient  response  of  combat  aircraft  skin-rib-strlnger  structures 
when  subjected  to  Internal  air  blast  loadings.  The  finite  elements  used  are 
flat,  isotropic,  rectangular  plates  and  beam  stiffeners.  The  input  pressure  Is 
assumed  constant  over  each  element.  The  output  consists  of  all  of  the  input 
data,  the  stable  time  increment,  the  stress-strain  data  for  the  panels  and  bars, 
and  the  deflections  at  the  node  points.  The  stresses  and  strains  are  tensor 
quantities  at  the  cross-section  at  the  center  of  the  elesnent.  Two  versions  of 
the  code  were  developed;  one  for  IBM  computers  and  one  for  CDC  computers.  The 
theory,  user's  manual  and  code  listing  of  the  CDC  version  are  given  in  Ref.  23. 
The  IBM  version  of  the  BR-1  code  was  modified  by  the  author  so  that  it  could 
be  used  to  predict  the  response  of  aircraft  fuel  tanks  when  subjected  to  the 
hydraulic  ram  pressure  loading  due  to  ballistic  penetrators.  The  modification 
was  based  upon  the  piston  theory  described  in  Section  III.  The  modified  code 
is  called  BR-lHR.  A description  of  che  theory,  the  modifications,  and  the 
additional  instructions  required  to  operate  BR-lHR  are  given  in  Ref.  24.  The 
code  is  operational  on  the  IBM  360/67  in  FORTRAM  IV,  Level  H. 

The  BR-lHR  code  can  he  uped  to  compute  the  response  of  the  NUC  exit  wall 
to  the  hydraulic  ram  pressure  loading.  One  finite  element  model  selected  to 
represent  the  20  inch  x 20  inch  exit  wall  rhown  in  Fig.  2 is  Illustrated  in 
Fig.  6.  This  model  has  49  elements  and  64  nodes.  The  reasons  for  selecting 
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this  model  were: 

1.  The  number  of  nodes  must  be  small  in  order  to  keep  the  execution 
time  to  a reasonable  figure. 

2.  The  strain  gage  locations  should  be  close  to  the  center  of  the 
element  since  thats  where  the  strains  are  printed. 

3.  The  edge  and  comer  elements  should  be  small  since  considerable 
edge  effects  are  anticipated. 

4.  This  model  can  be  further  refined  by  subdi\lding  the  large  element 
that  contains  the  exit  point  into  smaller  elements  in  order  to  allow 
a more  localized  pressure  loading. 

The  accuracy  of  this  model  for  linear,  elastic  response  is  evaluated  in 
the  appendix,  where  the  exact  solution  and  the  BR-lHR  solution  are  compared  for 
a simply  supported  plate  subjected  to  a unltorm  step  pressure.  The  response 
is  obtained  for  both  free  and  fluid  backed  plates.  The  essential  result  of  the 
accuracy  study  is  that  the  displacements  are  accurately  predicted  by  the  finite 
element  model  but  the  strains  are  not.  In  particular,  the  strains  at  the  center 
portion  of  the  plate  are  Inaccurate.  However,  it  should  be  noted  that  these 
strains  are  much  smaller  than  the  maximum  strains  in  the  plate.  The  solution 
for  a 144  node  model  was  also  obtained  in  order  to  determine  if  a smaller  p.rld 
size  would  yield  more  accurate  results.  The  144  node  model  provides  somewhat 
more  accurate  strains,  but  they  are  still  quite  inaccurate  in  the  center  portion 


of  the  plate. 
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The  BR-1  and  BR-IHR  codes  are  extremely  computer-time  consuming  codes*. 
Execution  times  on  the  IBM  360/67  for  small  to  moderate  size  problems 
(9-144  nodes)  vary  from  0.4  sec/node/time  step  to  0,6  sec/node/time  step, 
depending  to  some  extent  upon  the  amount  of  output.  Because  of  numerics' 
stability  requirements,  the  maximum  time  step  for  stable  results  is  very 
small  - of  the  order  of  microseconds. 

The  maximum  time  step  for  numerical  stabil Uy  for  the  64  node  model  of 
the  NWC  exit  wall  shown  in  Fig.  6 is  approximately  13,5  psec.  Furthermore, 
when  the  elements  of  the  model  are  made  smaller  for  better  solution  accuracy 
the  maximum  stable  time  step  decreases.  Thus,  due  to  the  very  small  allowable 
time  step,  the  large  number  of  nodes  required  for  accuracy,  and  the  large 
amount  of  time  required  per  node  at  each  time  step,  the  total  execution  time 
can  be  considerable.  For  example,  the  total  execution  time  required  to  com- 
pute the  response  of  the  64  node  model  for  A msec  at  12  psec/time  step  is  esti- 
mated to  be  approximately  10,000  secs  or  about  3 hours  on  the  IBM  3o0/67. 


*The  modifications  made  to  BR-1  did  not  significantly  effect  the  execution  time 
required  for  a solution.  Computer  execution  time  for  a problem  run  using  BR-IHR 
is  essentially  the  same  as  the  execution  time  using  BR-1. 
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B,  RATMS 


SATANS  (^tatlc  And  Transient  ^alysis,  Nonlinear,  ^hells)  is  a digital 
computer  code  that  computes  the  geometrically  nonlinear  static  and  dynamic 
response  of  arbitrarily  loaded,  elastic  shells  of  revolution.  SATANS  has 
been  in  existence  since  1967,  and  has  undergone  several  significant  modifica- 
tions since  the  static  version  was  first  described  in  Ref.  31.  The  latest 
public  version,  available  through  COSMIC  (LAR-11109)  can  be  used  to  analyze 
any  shell  of  revolution  for  which  the  following  conditions  hold: 

(1)  The  geometric  and  material  properties  of  the  shell  are 
axlsynaetrlc,  hut  may  vary  along  the  shell  meridian. 

(2)  The  shell  material  is  isotropic,  but  the  modulus  of 
elasticity  E may  vary  through  the  thickness. 

(3)  The  applied  pressure  and  temperature  distributions  can  be 
expressed  in  a Fourier  sine  or  cosine  series  in  the  cir- 
cumferential coordinate. 

(4)  The  boundaries  of  the  shell  may  be  closed,  free,  fixed  or 
elfistlcally  restrained. 

The  governing  partial  differential  equations  are  based  upon  Sanders’ 
nonlinear  thin  shall  theory  for  the  condition  of  f?mall  strains  and  moderately 
small  rotations.  The  Inplane  and  normal  inertial  forces  are  accounted  for, 
but  the  rotary  inertial  terms  are  neglected.  The  set  of  governing  nonlinear 
partial  differential  equations  is  reduced  to  an  Infinite  number  of  sets  of 
four  second-order  partial  differential  equations  in  the  meridional  and  time 
coordinates  by  expanding  all  dependent  variables  in  a sine  or  cosine  series 
in  terms  of  the  circumferential  coordinate  G . The  four  state  variables  are 
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the  three  displacements  and  the  meridional  bending  moment.  The  sets  are  un- 
coupled by  utilizing  appropriate  trigonometric  Identities  and  by  treating  the 
nonlinear  coupling  terms  as  pseudo  loads.  The  meridional  derivatives  are 
replaced  by  the  conventional  central  finite  difference  approximations,  and  the 
displacement  accelerations  are  approximated  by  the  implicit  Houbolt  backward 
differencing  scheme.  This  leads  to  sets  of  algebraic  equations  in  terms  of 
the  dependent  variables  z and  the  Fourier  Index  n.  Each  set  is  tridiagonal 
In  the  matrix  sense  with  4x4  matrices  on  each  of  the  diagonals.  Each  set 
of  algebraic  equations  is  solved  using  Potters*  form  of  Gaussian  elimination. 
For  considerably  more  detail  concerning  SATAMS  and  its  use  the  Interested 
reader  Is  referred  to  Refs.  25  and  26. 

In  SATANS,  the  applied  pressure  is  defined  at  each  finite  difference  sta- 
tion for  each  circumferential  harmonic,  and  the  output  consists  of  the  stress 
resultants  and  displacements  at  each  station  for  earh  circumferential  harmonic 
at  each  time  step.  For  this  study  SATANS  was  modified  to  accoun  • for  the 
piston  theory  radiation  damping  pew  in  the  transverse  equation  of  motion. 

The  output  was  also  changed  to  print  the  strains  at  each  station. 

SATANS  is  a very  fast  code  compared  to  BR-lHR.  For  example,  an  analysis 
of  the  axisymmetric  response  of  a flat  circular  plate  requires  approximately 
0.01  sec/station/ time  step.  Thus,  in  this  instance  its  approximately  50  times 
faster  than  BR-lHR*.  However,  it  cannot  account  for  inelastic  effects,  as 
BR-IHR  can. 

*The  axisytnnetric  analysis  uses  only  the  n ■ 0 mode.  An  asymmetric  analysis 
requires  additional  modes  and  hence  more  computational  time.  A reasonable 
comparison  for  asymmetric  resoonse  would  be  a 60  node  model  for  BR-lHR,  re- 
quiring 30  sec/time  step,  verses  a 15  station,  6 mode  model  for  SATANS, 
requiring  .9  sec/time  step. 


Since  the  NFS  entry  wall  is  a flat  circular  plate  with  a central  hole 
and  is  subjected  to  an  essentially  axisyoffletrlc  hydraulic  ram  pressure  loading, 
SATANS  provides  the  appropriate  analysis  for  purposes  of  comparison,  provided 
the  stresses  in  the  plate  remain  elastic.  However,  for  the  exit  wall  study 
SATAHS  provides  only  an  estimate  of  response,  since  the  plate  is  square,  the 
loading  is  not  symmetric,  and  the  plate  stresses  exceed  the  elastic  limit  near 
the  exit  point.  Nevertheless,  due  to  ':he  considerable  savings  in  computer  exe- 
cution time,  SATANS  is  used  for  purposes  of  comparison.  A typical  21  station 
model  of  the  exit  wall  is  shown  in  Fig,  6,  where  the  plate  has  been  centered  at 
the  exit  point  of  shot  14. 
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C.  Cavitation 

If  cavitation  takes  place  In  the  fluid  near  the  wall»  negative  Interface 
pressures  cannot  exist,  and  the  net  pressure  q drops  to  zero.  SATANS  has 
been  further  modified  to  account  for  cavitation  by  computing  the  net  Interface 


pressure  2p^  - pew  at  each  station,  at  each  time  step,  and  at  each  Iteration. 
If  the  net  pressure  at  any  station  Is  found  to  he  less  than  zero,  then  the  net 
pressure  at  that  station  la  set  to  zero  and  another  Iteration  Is  taken.  This 
In  effect  separates  the  wall  from  the  fluid  at  those  stations  where  cavitation 
occurs.  If  the  net  pressure  later  returns  to  positive  values  then  the  positive 
pressure  Is  applied  to  the  v;all.  This  procedure  takes  advantage  of  the  Itera- 
tive solution  procedure  used  by  SATANS  since  the  net  pressure  depends  upon  the 
velocity,  which  In  turn  depends  upon  the  pressure.  It's  not  possible  to  In- 
corporate cavitation  in  BR-IHR  exactly  since  an  explicit  time  integration  method 
Is  used  without  Iteration. 


V HYDRAULIC  RAM  PRESSURE  PREDICTIONS 
AND  COMPARISONS 

A.  The  Nl^C  HYDRAULIC  RAM  PRESSURE  CODE 

The  NWC  hydraulic  ram  pressure  code  Is  a digital  computer  code  developed 
to  predict  the  pressure  and  velocity  ttiroughout  a rectangular  volume  c£  fluid 
due  to  a ballistic  penetrator.  The  theory  of  the  code  and  Instructions  for 
its  use  are  deucrlted  in  detail  in  Ref.  5.  This  code  provides  toe  values  of 
the  hydraulic  ram  pressure  to  be  used  in  the  structural  response  codes  as  the 
incident  pressure  at  the  tank  wall. 

A simplified  model  of  the  complex  bullet  behavior  is  used  for  the  pressure 
computation.  The  bullet  is  assumed  to  travel  in  a straight  line.  The  rate  of 
kinetic  energy  loss  of  the  bullet  is  taken  to  be 

■ m0V^ 

where  E is  the  kinetic  energy  of  the  bullet,  Xj^  is  the  bullet  position, 
m is  the  bullet  mass,  V is  the  bullet  velocity,  and  S is  a velocity  decay 
coefficient  that  is  a function  of  x^  for  tumbling  bullets.  In  the  code  the 
bullet  is  assumed  to  enter  the  test  cell  with  0-degree  yaw  and  to  continue  in 
this  attitude  until  the  distance  XTUM  is  reached  where  it  begins  to  tumble. 

The  bullet  is  assumed  to  be  fully  tumbled  at  a distance  XTUM  plus  DXTUM 
Stripping  of  bullet  jackets  from  the  atmor  piercing  core  of  API  ammunition  is 
also  accounted  for,  albeit  in  a crude  manner. 

The  wave  equation  is  used  to  calculate  the  flow  field  resulting  from  the 
bullet  and  cavity  motion.  The  projectile  and  cavity  are  approximated  by  the 
action  of  a line  of  sources  along  the  bullet  path.  The  stress  waves  are 
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assumed  to  propagate  in  a spherical  manner  from  each  source.  The  source 
strength  distribution  is  estimated  by  equating  the  energy  loss  of  the  bullet 
to  the  sum  of  the  work  done  by  the  difference  in  ambient  and  cavity  pressure 
and  the  kinetic  energy  of  the  fluid.  The  pressure  in  the  fluid  is  computed 
using  Bernoulli’s  equation,  Eq.  (2).  Although  the  action  of  the  cavity 
Is  accounted  for  in  the  fluid  model,  the  absence  of  fluid  within  the  cavity 
is  ignored.  Thus,  the  code  pressures  for  locations  within  the  cavity  are  of 
dubious  value.  This  aspect  is  discussed  in  the  following  section. 

An  important  feature  of  the  code  is  its  ability  to  account  for  reflec- 
tions from  each  of  the  six  surfaces  of  the  fluid  volvime  using  the  method  of 
Images.  Two  options  are  available;  either  reflections  from  free  surfaces,  or 
reflections  from  rigid  walla*.  Since  the  piston  theory  requires  either  twice 
the  incident  hydraulic  ram  pressure  (2p^)  or  the  pressure  at  a rigid  wall  ( ^e 
two  being  the  same  when  the  velocity  squared  term  is  neglected  in  Bernoulli’s 
equation),  either  option  can  be  used  to  obtain  the  pressure  at  the  wall.  However, 
since  the  other  walls  of  the  tank  are  not  rigid;,  only  reflections  from  the 
wall  being  considered  should  be  accounted  for  if  the  rigid  wall  approach  is 
used**.  Examples  aud  users  instructions  on  the  image  feature  are  given  in 
Ref.  5. 


*The  free  surface  boundary  condition  p « 0 is  satisfied  for  the  potential  term 
only.  In  order  to  have  the  rigid  wall  refelctions,  the  statement  ASSIGN 
(I)  « 1.0  must  be  inserted  after  530  in  SUBROUTINE  IMAGE. 

**Mo8t  aircraft  fuel  tank  walls  probably  reflect  more  like  free  surfaces  than 
rigid  walls.  Refs.  3 and  18. 
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The  input  parameters  to  the  code  include  the  dimensions  of  the  rectangu- 
lar volume  of  fluid  into  which  the  penetrator  is  fired,  XC(l)  - XC(3) , and  the 
entrance  and  exit  points  of  the  projectile,  X(l)  - X(6) . in  order  to  obtain 
the  hydraulic  ram  pressure  at  the  wall  being  considered  the  user  must  specify 
the  locations  of  the  points  at  which  these  variables  are  to  be  computed  and  the 
time  increment  for  which  the  pressure  is  to  be  computed.  Thus,  in  the  BR-IHR 
code  the  center  of  each  element  is  the  location  for  the  pressure  computatirn, 
whereas  in  SATANS  the  location  of  the  finite  difference  stations  defines  the 


points  for  the  pressure  computation*. 

Additional  input  data  to  the  NWC  code  is  given  below: 

XMASS(l)  - bullet  mass,  lb. 

XMASS(2)  - penetrator  mass,  lb. 

XMASS(3)  - Jacket  mass,  lb. 

2 

AR£A(1)  - bullet  presented  area,  normal  attitude,  in 

2 

AR£A(2)  - bullet  presented  area  tumbled  attitude,  in 

2 

AREA(3)  - penetrator  presented  area,  normal  attitude,  in 

2 

AB£A(4)  - penetrator  presented  area,  tumbled  attitude,  in 

2 

AR£A(5)  - Jacket  presented  area,  tumbled  attitude,  in 

2 

AREA(6)  - bullet  stern-first  presented  area,  in 
DRAG(l)  - coefficient  of  drag  for  bullet,  normal  attitude 

DRAG(2)  - coefficient  of  drag  for  bullet,  tumbled  attitude 

DRAGC3)  - coefficient  of  drag  for  penetrator,  normal  attitude 

DRAG(A)  - coefficient  of  drag  for  penetrator,  tumbled  attitude 


*The  hydraulic  ram  pressure  is  either  2p  or  the  reflected  pressure  at  a rigid 
wall.  Since  the  NWC  code  does  not  proviSe  an  option  to  punch  the  pressure,  the 
code  was  modified  by  the  author  to  punch  these  data  in  the  format  required  by 
BR-lHR.  This  required  changes  to  the  subroutine  MIRROR  in  the  procedure  for 
computing  a minimum  initial  time.  These  changes  established  a common  minimum 
initial  time  to  be  used  at  all  points.  The  punch  routine  provides  the  pressure 
the  wall  at  each  time  increment,  starting  with  the  minimum  initial  time. 
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DRAG (5) 

DRAG(6) 

DENS 

PO 

PC 

c 

BC 

VKL 

DX 

NT 

XSTRIP 

XTUM 

nxTim 

NE 

XE 

XMAX 

EN 

EFRACT 

NCTUM 

BS,0MEG 

NSTRIP 


coefficient  of  draj?  for  jacket,  tumbled  attitude 
coefficient  of  drag  for  bullet  stern-first,  normal  attitude 
fluid  density,  Ib/ln^ 
ambient  pressure,  psi 

- cavity  pressure,  psi,  normally  taken  to  be  zero 

sound  speed  in  fluid,  ft/sec 

constant,  normally  taken  to  be  0.434 

initial  bullet  velocity,  ft/sec 

distance  Increment  along  trajectory,  in, 

a variable  such  that  NT  - 0/1  means  that  the  bullet  does 
not  tumble/does  tumble 

distance  along  the  trajectory  that  the  bullet  begins  to  strip 

distance  along  the  trajectory  that  the  bullet  begins  to  tumble, 
in* 

- distance  required  for  the  bullet  to  become  fully  tumbled,  in, 
if  NT  » 0,  xnjM  and  DXTUM  can  be  arbitrary 

- is  a variable  such  that  NE  - 0/1  means  bullet  does  not  exit 
tank/does  exit  tank 

- distance  along  the  trajectory  where  the  bullet  exits  the 
volume,  in.  If  NE  ■ 0,  XE  can  be  arbitrary 

- maximum  length  of  trajectory  through  volume,  in.  Dimension 
statements  limit  XMAX  to  300*DX. 

- exponent  which  specifies  power  law  followed  by  the  drag 
function  (EN  ■ 3.0  in  general) 

- a factor  related  to  bullet  stripping  (EFRACT  « 3.0  in  general) 

- a variable  such  that  NCTUM  •»  0,1,2  - means  bullet  does  not 
continuously  tumble/bullet  continuously  tumbles/bullet  con- 
tinuously tumbles  for  one  cycle 

- parameters  not  used  for  this  version,  set  equal  to  1 

- a variable  such  that  NSTRIP  “0,1  means  bullet  does  not 
strip/ballet  does  strip 
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The  code  Is  programmed  such  that  the  computation  of  the  trajectory  para-' 

meter  Is  automatically  stopped  when  the  bullet  reaches  the  exit  wall  (XMAX) . 

However,  its  possible  to  compute  pressures  for  times  longer  than  the  time  the 

bullet  takes  to  penetrate  the  tank.  Since  the  theoretical  model  asstunes 

spherical  propagation  from  a line  of  sources,  the  source  strength  theoretically 

drops  to  zero  when  the  bullet  exits  the  tank.  The  program  indicates  this  fact, 

when  no  reflections  are  considered,  by  printing  a constant  value  for  (PPHI) 

at  all  times  .'.ater  than  the  time  required  for  the  zero  source  disturbance  to 

reach  the  point  in  question.  However,  the  pressure  (P)  printed  at  these  times 

2 

(Eq.  2)  cnanges  due  to  the  change  in  V , which  is  based  upon  the  invalid 
trajectory  parameters.  When  reflections  are  considered,  <f>  also  changes  due 
to  the  image  sources. 

B.  The  Cavity 

In  the  series  of  ballistic  tests  conducted  by  NWC  to  obtain  the  detailed 
pressure  measurements  (Ref.  4),  highspeed  motion  pictures  were  taken  through  a 
one  inch  plexiglass  window  at  the  side  of  the  tank.  This  provided  a 30  inch 
high  by  36  inch  long  field  of  view*.  The  opposite  plexiglass  wall  was  sand- 
blasted and  thus  acted  as  a diffusing  screen  for  backlighted  photography.  Th® 
films  vividly  show  the  bullet  penetrating  the  entry  wall  and  passing  through 
the  fluid,  the  resulting  cavity  caused  by  the  bullet,  and  the  entry  wall  deflec- 
tion. Fig.  7 presents  several  frames  of  the  film  showing  the  penetration  by  a 
30  caliber  projectile**.  Several  features  regarding  the  cavity  can  be  seen  from 
these  pictures.  First,  the  cavity  does  not  touch  the  entry  wall  anywhere  except 
in  the  immediate  vicinity  of  the  entry  hole.  Second,  there  does  not  appear  to 


*The  tank  was  a five  foot  cube. 

**The  pipes  extending  into  the  tank  from  the  top  contain  the  pressure  trans- 
ducers. The  first  pipe  is  6 inches  from  the  entry  wall  and  the  pipes  are 

spaced  6 inches  apart. 
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b«  any  cavitation  In  the  fluid  at  the  entry  wall  due  to  negative  interface 
pressures  caused  either  by  the  wall  deflection  (away  from  the  fluid)  or  by 
tension  waves  In  the  fluid.  However,  pictures  of  later  tines  definitely  show 
cavitation  on  the  side  walls.  Third,  when  the  bullet  has  traveled  approximately 
24  Inches  Into  the  tank  the  cavity  Is  conical  for  about  12  Inches,  followed  by 
a larger,  more  cylindrical  cavity.  The  diameter  of  the  cylindrical  part  of  the 
cavity  Is  approximately  six  Inches.  When  the  bullet  has  progressed  approximately 
30  inches,  the  cylindrical  cavity  haa  grown  to  about  eight  Inches.  The  size  of 
the  cavity  Is  Important  when  attempting  to  determine  the  net  pressure  on  the  exit 
wall  after  the  bullet  exits  the  tank,  l.e.  vdiere  Is  there  contact  between  the 
fluid  and  the  wall? 
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C.  Entry  Wall  Pressure 

The  loading  on  the  entry  wall  of  the  fuel  tank  due  to  the  ballistic  peni — 
trator  can  be  divided  into  four  successive  phases: 

1.  a penetration  phase  where  the  projectile  punches  a hole  in  the 
entry  wall, 

2*  a shock  phase  where  the  projectile  protuslon  into  the  fluid  causes 
a hemispherical  shock  wave, 

3.  a drag  phase  during  which  the  kinetic  energy  of  the  projectile  is 
dissipated  by  form  drag  and  a cavity  behind  the  projectile  is  formed, 

4.  a cavity  oscillation  phase  which  results  from  the  growth  and  collapse 
of  Idle  cavity. 

The  penetration  phase  was  not  Included  in  the  NFS  tests  since  the  entry  plate 
had  a one  inch  diameter  hole  at  the  center  through  which  the  projectile  passed. 
This  was  done  to  eliminate  the  first  phase  from  the  problem. 

The  NFS  prediction  of  the  fluid  pressure  during  the  shock  phase  uses  the 
Yurkovich  theory  which  is  based  upon  the  assumptions  of  a point  energy  source 
release,  a rigid  entry  wall,  and  a shock  radius  that  is  proportional  to  a con~ 
stant  power  of  time.  A typical  predicted  wave  front  pressure  distribution 
resulting  from  the  impact  of  a 45  grain  .222  caliber  bullet  into  a water  filled 
tank  is  shown  in  Fig.  8.  Note  thr.t  this  phase  is  much  less  than  one  msec  in 
duration.  The  magnitude  of  this  phase  relative  to  the  drag  phase  is  heavily 
dependant  upon  the  relative  amount  of  kinetic  energy  lost  by  the  projectile 
during  the  shock  phase,  which  in  turn  is  heavily  dependant  upon  the  projectile 
Mach  number  with  respect  to  the  fluid.  Projectiles  with  high  Mach  numbers  (>2) 
lose  considerably  more  energy  in  the  shock  phase  than  those  wltVi  low  Mach 
numbers.  Detailed  shock  phase  results  from  the  NFS  experiments  and  comparisons 
with  the  results  from  the  Yurkovich  theory  are  given  in  RefLi.  9,  12,  and  13. 
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Since  the  shock  nhase  does  not  dissipate  all  of  the  available  projectile 
kinetic  energy,  the  projectile  travels  through  the  tank  until  either  all  of  the 
energy  is  dissipated  or  it  strikes  and/or  penetrates  the  exit  wall.  No  pro- 
jectiles penetrated  the  exit  wall  in  the  NPS  tests.  The  predicted  fluid 
pressures  on  the  entry  wall  during  the  drag  phase  can  be  obtained  from  the 
NWC  hydraulic  ram  pressure  code. 

Values  of  the  input  data  used  in  the  NWC  hydraulic  ram  pressure  code  are  given 
below: 


EN 

BS 

EFRACT 

NCTUM 

NS TRIP 

NT 

XTUM 

DXTUM 

XSTRIP 

XMASS(3) 

AR£A(6) 

DRAG (6) 

BETA(7) 

VEL 

DENS 

PO 

PC 

C 

BC 


3.00000 

1.00000 

3.00000 

0 

0 

0 

20.00000 

25.00000 

20.00000 

0.00640  0,00010  0,00010 

0.03800  0.03800  0.03800  0.03800  0,03800  0.03800 

3.00000  3.00000  3.00000  3,00000  3.00000  3.00000 

0.32063  0.32063  20.51999  20.51999  20.51999  0.32063  20.51999 

2500.00000 
0.03600 

14.70000 

0.0 

4900.00000 
0.43400 


Note  that  tumbling  of  the  bullet  was  not  considered  since  no  data  on  tumbling  be- 
havior was  available.  Instead,  a constant  drag  coefficient  for  the  bullet  was 
determined  from  the  projectile  position  data.  Neglect  of  the  ttmbling  behavior 
strongly  influences  the  shape  of  the  pressure  curves.  Figure  9 shows  comparisons 
of  the  measured  and  predicted  fluid  pressure  2 in.,  3 in.,  ond  4 in.  from  the 
entry  point  along  a radius  50”  from  the  vertical  and  in  the  direction  of  the 
projectile  trajectory.  The  energy  level  was  7,493  in. -lb  and  the  entry  wall 
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was  0,5  in,  thick.  Rigid  wall  reflections  from  all  surfaces,  including  the 
front  wall,  were  assumed.  The  predicted  shock  pressures  are  included  in  the 
comparison.  The  shock  pressures  were  obtained  based  upon  the  assumption  of 
a 10»  loss  of  projectile  energy  during  the  shock  phase.  The  measured  and  pre- 
dicted fluid  pressure  at  2 in. , 4 in.  and  6 in.  for  the  condition  of  free 
surface  every\^ere  is  shown  in  Fig.  10.  The  shock  phase  pressures  are  not  in- 
cluded in  this  figure.  Examination  of  Figs.  9 and  10  reveals  that  the  predicted 
drag  phase  pressure  at  the  larger  radii  have  essentially  the  same  peak  pressure 
as  the  measured  pressure  but  the  wave  form  is  different.  The  difference  in  form 
may  be  due  to  the  fact  that  the  tumbling  behavior  is  not  properly  accounted  for. 
Furthermore,  note  that  the  measured  pressure  trace  does  not  Indicate  a shock 
phase.  Thl';  may  be  due  to  the  limitations  on  the  time  rise  of  the  charge  ampli- 
fier used*.  Much  more  detailed  Informtlon  concerning  the  fluid  pressure 
measurements  near  the  entry  wall  is  given  in  Refs.  8,  9,  11,  12,  13,  14,  16, 
and  17. 

The  predicted  hydraulic  ram  pressure  on  the  entry  wall  for  the  7,493  in. -lb 
energy  level  shot  is  shown  in  Fig.  11  at  SOysec,  lOOysec,  and  ISOysec  after 
impact..  The  hydraulic  ram  pressure  is  taken  as  the  sum  of  twice  the  incident 
drag  pressure  and  the  shock  pressure.  Free  surface  reflections  from  every 
surface  except  the  entry  wall  were  included**.  This  transient  pressure  field 
is  the  one  selected  for  use  with  SATANS. 


*Shadovgraphs  show  the  existence  of  a hemispherical  shock  front.  Ref.  13. 

**The  choice  was  made  to  use  twice  the  incident  drag  phase  pressure  rather  than 
the  reflected  pressure  at  a rigid  wall  since  this  allows  the  use  of  free  surface 
at  the  other  walls.  A comparison  of  the  two  pressures  on  one  shot  revealed  a 
difference  of  less  than  4%. 
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Prcrssura  (psi)  Pressure  (psi)  Pressure  (psi) 


XXX  Experimental 
Calculated  Drag  Phase 


Rodius  « 2ln 


Time  (ft sec) 

FIG.  IQ  FLUID  PRESSURE  COMPARISON,  FREE  SURFACES 
(REF  17) 


D.  Exit  Wall  Pressure 

The  hydraulic  ran  pressure  on  the  exit  wall  can  also  be  predicted  using 
the  NWC  code. 

Values  of  the  Input  data  for  shot  14  of  the  NWC  test  program  are  given  below: 


EN 

3.00000 

BS 

i.oqooo 

EFRACT 

3.00600 

NCTUM 

1 

NSTRIP 

0 

NT 

1 

XTUM 

0.67929 

DXTUM 

11.35918 

» 

XSTRIP 

,35.00000 

XMASSC3) 

0.10660 

0.06410 

0.04250 

AREA(6) 

0,20460 

1.03780 

0.14320 

0.70020 

1.03780 

0.08550 

DRAG(6) 

0.05000 

0.30000 

0,05000 

' 0.30000 

1.00000 

0.82000 

BETA(7) 

0.00173 

0.05272 

0.00202 

0.05915 

0.44076 

0.01187 

BETA(7)  - 0.00173  0.05272  0.00202  0.05915  0.44076  0.01187  C. 01974 

The  tumbling  parameters  XTTJM  and  DXTUM  were  obtained  by  determining  the  best  fit 
of  the  predicted  pressure  to  the  measured  pressure  at  one  point  In  the  fluid. 

The  entry  point  of  the  projectile  was  30  In.,  30  In.,  and  0 in.  In  the  x,  y,  and 
z directions  respectively.  The  exit  point  was  28  In.,  2.5  In.  and  22  In.  The 
projectile  was  predicted  by  the  code  to  strike  the  exit  wall  at  1.63  msec  after 
entering  the  tank.  The  pressure  build  up  started  at  0.35  msec.  This  tine  is 
taken  as  the  zero  time  for  the  structural  response  computations . 

Two  different  Image  conditions  are  considered  for  shot  14  In  order  to  Illus- 
trate the  effect  of  reflections  from  the  other  walls  of  the  tank;  that  of  no 
image  and  that  of  free  surface  reflections  from  every  wall  except  the  exit 
wall.  The  predicted  pressure  2p^  at  the  center  of  the  exit  wall,  at 
X > y ■ 26.5  in.,  and  at  the  exit  point  Is  given  In  Fig.  12  as  a function  of 
time  for  the  condition  of  no  Images  and  In  Fig.  13  for  the  condition  with  free 
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Time  (msec) 


FIG.  12  EXIT  WALL  PRESSURE  2p.,N0  IMAGES. SHOT  14 


pressure  drops  tozi 


(psi) 


Piai3  EXIT  WALL  PRESSURE  2g,  FREE  SURFACE  REFLECTIONS 


Pressure  drops  to  zero 


surface  reflections.  The  pressures  given  In  Figs.  12  and  13  ere  assumed  to 
drop  to  zero  shortly  after  the  bullet  exits  the  tank,  for  the  reasons  given 
In  Section  V.A,  The  pressure  field  2p^  on  the  exit  wall  at  t “ 0.9  msec, 

1.6  msec,  and  1.75  msec  Is  shown  In  Figs.  14-16  for  the  conditions  of  no 
Images,  and  In  Figs.  17-19  for  the  condition  of  free  surface  Images.  The 
region  around  the  exit  point  where  the  pressure  Is  assumed  to  be  zero.  The 
zero  pressure  there  Is  due  to  the  fact  that  the  energy  source  has  exited  the 
tank;  It  Is  not  related  to  the  cavity  that  follows  the  bullet.  However,  the 
presence  of  the  cavity  will  also  cause  the  pressure  to  drop  to  zero. 

f 

A comparison  of  the  pi;-easures  on  the  exit  wall  without  reflections  with 
the  pressures  with  reflections  reveals  a considerable  difference  between  the 
two  pressure  fields,  particularly  at  times  later  than  1.5  msec.  In  general, 
the  no  Image  pressures  are  higher  than  those  with  Images.  Figure  18  shows 
negative  pressures  over  a considerable  portion  of  the  plate  at  1.6  msec.  In 
addition.  Fig.  13  shows  some  strong  pressure  spikes  Just  after  the  bullet  exits 
the  tank.  Its  not  clear  where  these  spikes  come  from. 

As  a consequence  of  both  the  doubtful  validity  of  the  predicted  fluid 
pressure  on  the  wall  after  the  bullet  exits  the  tank  and  the  presence  of  the 
cavity,  the  assumption  Is  made  in  the  analysis  for  the  exit  wall  response  that 
the  pressure  field  2p^  drops  to  zero  when  the  bullet  exits  the  tank.  This 
assumption  of  course  neglects  the  Incoming  pressures  that  have  been  reflected 
from  the  other  surfaces,  but  the  presence  of  the  cavity  makes  it  questionable 
whether  or  not  the  fluid  Is  In  contact  with  the  wall.  High-speed  films  of  the 
exit  wall  show  that  the  water  does  not  Immediately  pour  out  of  the  exit  hole 
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FIG.  17  EXIT  WALL  PRESSURE  2^,  FREE  SURFACE 
REFLECTIONS,  t=  0.9msec 
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after  the  bullet  penetrates  the  wall.  There  appears  to  be  some  delay,  indicating 
the  1 esence  of  the  cavity  In  the  vicinity  of  the  hole.  Later,  however,  the 
water  is  periodically  pumped  out  in  large  quantities  during  the  oscillating 
cavity  phase. 
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VI  WALL  STRAIN  COMPARISONS 


A,  Entry  Wall  Strains 

Figure  20  shows  typical  experimental  entry  wall  strains  for  a 0,16 
Inch  thick  aluminum  wall.  Additional  data  are  contained  In  Refs.  IS  and  17. 

The  upper  photograph  shows  that  the  radial  strain  Is  essentially  due  to  pute 
bending  because  of  the  relatively  symmetrical  traces.  The  lower  photograph 
shows  that  circumferential  strains  are  also  primarily  due  to  pure  bending  of 
the  plate  for  the  first  400  iiseconda.  However*  after  400  pseconds  the  circum- 
ferential strains  are  not.  as .symmetric*  which  Is  Indicative  of  stretching  of 
the  middle  surface  of  the  plate.  The  traces  for  the  0.09  In.  plate  are  similar 
to  tbps^  shown  In  Fig.  20*  but  they  reveal  more  stretching  of  the  middle 
surface. 

Strains  in  the  0.09  in.  entry  wall  have  been  computed  using  SATANS,  The 
applied  transient  pressure  field  on  the  wall  was  the  shock  and  cavity 
pressure  shown  In  Fig.  11,  The  7.5  In.  radius  plate  was  modeled  with  71 
equally  spaced  nodes,  starting  at  the  inner  radius  of  0,5  In,  and  continuing 
outward,  0.1  In.  apart.  The  boundary  conditions  at  the  inner  hole  were  those 
of  a free  edge;  fully  clamped  conditions  were  assumed  at  the  outer  edge.  The 
time  increment  was  5 ysec.  Poisson's  ratio  was  taken  as  0.33. 

Several  cases  of  fluid-structure  interaction  were  considered.  When  the 
effect  of  the  fluid  on  the  wall  motion*  i.e.  the  -pc<>  term  in  Eq.  (7),  was 
ignored,  the  computed  strains  caused  by  the  sum  of  the  shock  and  drag  pressures 
were  very  large  compared  to  the  measured  strains.  When  the  piston  theory  Is 
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used,  the  predicted  strains  are  reduced  by  a considerable  amount.  Figure  21 
presents  the  predicted  and  the  measured  radial  strains  on  both  sides  of  the  wall 
at  2 In.  from  the  entry  point  as  a function  of  time.  Figure  22  shows  the 
corresponding  circumferential  strains* **. 

The  sharp  rise  and  fall  of  the  predicted  radial  strain  around  50  ysec  in 
Fig,  21  is  probably  due  to  the  passage  of  the  shock  pressure  front.  Note  that  the 
the  early  strains  Indicate  compression  on  the  dry  side  and  tension  on  the  wet 
side.  The  reason  for  this  can  be  seen  by  examining  a plot  of  the  predicted 
wall  deflection  at  SO  ysec  shown  in  Fig.  23. 

In  general,  the  form  of  the  predicted  strains  In  Figs.  21  and  22  Is  the 
same  as  the  form  of  the  measured  strains.  However,  except  for  the  radial 
strains  due  to  the  shock  pressure,  the  predicted  strains  are  considerably 
Smaller  than  the  measured  strains.  One  possible  exv>lanatlon  for  this  difference 
Is  that  cavitation  takes  place  In  regions  where  the  net  Interface  pressure  is 
negative  (tension)*.  The  predicted  net  Interface  pressure,  2p^  (drag)  4p^ 
(shock)  -pew  , Is  plotted  In  Fig.  24  at  50,100  and  150  psec  after  Impact.  This 
figure  can  be  compared  with  Fig.  11  to  illustrate  the  effect  of  the  wall  motion 
on  the  Interface  pressure.  Note  the  presence  of  large  negative  pressures  in 
Fig.  24  The  response  of  the  entry  wall  assuming  cavitation  takes  place  has 
been  computed  using  the  modified  version  of  SATANS  discussed  In  Section  IV.  C. 


*A  similar  comparison  of  strains  is  presented  In  Ref.  17,  However,  In  Ref.  17 
the  sign  of  the  strain  appears  to  be  reversed  from  normal  convention,  l.e. 
tension  strain  Is  called  negative  and  compression  strain  Is  called  positive. 
This  conclusion  Is  based  upon  the  assumption  that  In  general  the  plate  will 
stretch,  not  compress,  after  the  initial  wavefront  passes.  In  addition,  the 
shock  pressure  used  in  Ref.  17  was  multiplied  by  a factor  of  two. 

**Recall  that  the  photographs  of  the  entry  plate  presented  in  Fig.  7 do  not 
appear  to  indicate  any  cavitation. 
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(b)  Dry  Side 


Strains  2 in.  From  Center  of  Plate 


o 


Predicted 

Measured 


(b)  Dry  Side 


Fig.  22.  Entry  Wall  Circumferential  Strains  2 in.  Frcm  Center  of  Plate 


R^essurft,  psi  Pressurt,  psi 


The  predicted  and  measured  results  for  the  radial  and  circumferential  strains 
on  both  sides  of  the  wall  at  2 in.  from  the  entry  point  are  presented  in 
Figs.  25  and  26.  Examination  of  these  results  reveals  that  the  strains  with 
cavitation  included  are  significantly  larger  than  those  computed  without 
cavitation.  However,  except  for  the  circumferential  strain  on  the  dry  side, 
the  general  form  of  the  predicted  strains  with  cavitation  differs  from  that 
of  the  measured  strains.  There  is  considerably  more  oscillation  in  the  pre- 
dicted strains.  This  is  because  the  damping  effect  of  the  fluid  has  been 
removed  in  those  areas  where  cavitation  occurs.  Thus,  perhaps  cavitation  is 
not  occurring  at  the  entry  wall.  Recall  that  the  NWC  high-speed  motion  pictures 
did  not  reveal  any  cavitation  at  an  entry  wall. 


(b)  Dry  Side 


Fig.  25.  Entry  Wall  Badial  Strains  2 in.  From  Center  of  Plate 
(with  Cavitation) 
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Predicted  / 

Measured  j 


(b)  Dry  Side 

Fig.  26.  Entry  Wall  Circumferential  Strains  2 in.  From  Center  of  Plate 
(with  Cavitation) 


H.  Exit  Wall  Strains 


As  an  example  of  the  strain  data  obtained  by  NWC  In  their  exit  wall  tests* 
the  strains  measured  at  the  16  gages  are  presented  for  shot  14  In  Fig.  27. 
Figure  28  shows  the  correspondence  between  the  gage  number  and  the  location  and 
direction  of  strain  measurement.  The  gage  locations  are  also  Indicated  In 
Figs.  2 and  6.  The  locations  defined  as  a and  b In  Fig.  27  at  gages  9^12  and 
13-16  respectively  are  also  Indicated  in  Figs.  6 and  28.  Note  in  Fig.  27 
that,  only  gage  6 appears  to  be  clearly  inoperative.  However,  the  results  from 
gage  16  are  also  suspect  since  they  are  of  a different  nature  than  those  from 
the  other  gages.  On  most  of  the  other  shots  many  more  gage  readings  were 
lost  because  the  wires  leading  from  the  gages  to  the  recording  device  were 
broken  by  the  exiting  bullet,  the  large  wall  motion,  and  the  spewing  water. 

As  a consequence  of  the  experimental  set-up,  the  time  coordinate  on  each 
strain  trace  has  a different  origin  with  respect  to  universal  time;  hence,  the 
strain  traces  cannot  be  compared  using  the  printed  time  as  a reference.  In 
addition,  there  Is  no  precise  Indication  or  data  on  the  exact  time  when  the 
bullet  strikes  the  exit  wall*.  Thus,  its  difficult  to  directly  compare  pre- 
dicted strains  with  the  measured  strains  on  a timewise  basis.  Nevertheless, 
some  indication  of  the  applicability  of  the  analysis  can  be  obtained  by  com- 
paring the  general  form  of  the  predicted  and  measured  strains  and  the  maximum 
strains . 


*An  examination  of  the  strain  at  gage  14  in  Fig.  27  reveals  a significant 
change  In  behavior  between  10.2  and  10.7  msec.  That  could  be  an  Indication 
that  the  bullet  exited  the  tank  at  that  time.  Furthermore,  a study  of  the 
high-speed  films  of  shot  12  revealed  that  two  of  the  wires  from  gages  9 
through  12  were  cut  by  the  bullet  about  four  frames  after  the  bullet  exited 
the  tank.  Thus,  the  signals  from  these  gages  were  probably  lost  approximately 
0.5  msec  after  penetration.  The  strain  trace  for  gage  9 was  definitely  lost 
at  about  5.7  msec.  Prior  to  that  time  the  strain  was  less  than  1000  ystrain. 
Since  the  maximum  strain  of  shot  12  was  approximately  4500  ystrain,  and  since 
gage  9 was  as  close  to  the  exit  point  as  any  of  the  other  gages.  It  appears 
that  the  large  strains  In  the  wall  occur  after  penetration,  except  of  course. 
In  the  immediate  vicinity  of  the  exit  point. 
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Fig.  27.,  Measured  Exit  Wall  Strains,  Shot  14  (Ref.  6) 


The  response  of  the  exit  wall  to  shot  14  has  been  computed  using  SATANS 

with  the  21  station  model  shown  in  Fig.  6.  The  stations  are  0.5  in.  apart. 

The  plate  is  assumed  to  be  continuous  at  the  center  (the  exit  point) , and  the 

outer  edge  is  assumed  to  be  fully  clamped.  The  pressure  2p  from  the  MC 

o 

hydraulic  ram  pressure  code  for  shot  14  along  the  line  x •=  28  in.  and 
centered  at  the  exit  point  was  applied  for  the  interval  t ■ 0.35  msec  to 
t • 1.625  msec  (fluid  time)  or  t •»  0 to  t “ 1.275  msec  (structure  time)  in 
increments  of  25  ysec*.  Tlie  pressure  was  assumed  to  be  zero  after  1.275  msec 
for  the  reasons  given  in  Section  V.D.  The  axisymmetric  response  of  the  circu- 
lar plate  was  computed  at  intervals  of  25psec  for  the  interval  t * 0 to  4 
msec.  The  results  for  the  radial  strain  (fy)  history  and  circumferential 
strain  (e^)  history  on  both  sides  of  the  wall  at  point  c in  Fig.  6 are  pre- 
sented in  Fig.  29  for  the  condition  of  no  pressure  reflections  from  the  other 
walls  and  in  Fig.  30  for  the  condition  of  free  surface  reflections.  A compari- 
son of  Figs.  29  and  30  reveals  that  the  strains  computed  with  the  reflections 
are  smaller  than  those  without  reflections.  Only  on  the  wet  side  shows 

any  major  change  in  the  form  of  strains. 

The  predicted  strain  traces  in  Figs.  29  and  30  can  reasonably  be  compared 
with  these  of  gages  9 through  12  (point  a in  Fig.  6)  or  13  through  16  (point 
b in  Fig.  6)  shown  in  Fig.  27.  The  time  scale  for  the  predicted  strains  has 
been  made  essentially  the  same  as  that  for  the  trt^asured  strains  for  ease  of 
comparison,  but  the  exact  times  cannot  be  compared.  The  location  of  each 
strain  in  Figs.  29  and  30  corresponds  to  the  strains  as  they  are  located  in 
Fig.  27. 

*The  minimum  time  at  which  the  fluid  pressure  starts  to  build  up  on  the  exit 
wall  is  0.35  msec.  Hiis  time  is  used  as  the  zero  time  for  the  structural 
response  computations. 
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A comparison  of  the  predicted  and  measured  strains  shows  that  the  shape 
ot  the  measured  traces  for  e on  the  dry  side  (gages  11  and  15) , e on 

X X 

the  wet  side  (gage  12;  gage  16  appears  to  be  Inoperative)  and  e^,  on  the  wet 
side  (gage  10)  is  very  similar  to  the  predicted  strain  history  for  those 
three  strains.  In  particular,  note  the  build  up  of  on  the  dry  side  and 

the  rapid  change  of  direction  of  e on  the  wet  side  when  the  bullet  exits 
the  tank.  However,  the  magnitude  of  the  predicted  strains  is  off  by  factors 
of  60,  200  and  30  respectively.  Furthermore,  the  predicted  strains  do  not 
show  the  high  frequency  content  exhibited  by  the  measured  strains.  This  could 
be  due  in  part  to  the  fact  that  the  analysis  is  axlsymmetric  and  does  not  con- 
tain the  high  frequency  asymmetric  motion  generated  by  the  corners  of  the 
square  plate,  and  in  part  to  the  fact  that  the  piston  theory  without  cavitation 
causes  the  plate  to  respond  as  a heavily  damped  oscillator.  An  examination 
of  the  net  fluid  pressure  on  the  exit  wall  revealed  that  negative  (tension) 
pressures  of  approximately  20  psi  or  less  occurred  over  portions  of  the  wall 
at  various  times.  To  Illustrate  this.  Fig.  31  presents  the  net  Interface 
pressure  over  the  wall  at  1,  2 and  3 msec.  When  the  arsumptlon  is  made  that 
cavitation  occurs  at  any  negative  pressure,  the  modified  version  of  SATANS 
predicts  the  strains  shown  in  Fig.  32  for  the  condition  of  no  reflections. 

A comparison  of  these  strains  with  the  measured  strains  (Fig.  27)  and  with  the 
predicted  strains  without  cavitation  (Fig.  29)  reveals  that  the  predicted 
strains  with  cavitation  are  much  closer  in  magnitude  to  the  measured  strains, 
but  the  similarity  in  form  has  been  lost. 

An  analysis  with  cavitation  using  a 30  in.  diameter  plate  centered  at 
the  exit  point  of  shot  14  gave  roughly  the  same  results  as  the  analysis  of 
the  20  in.  plate. 
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A general  comparison  of  maximum  strains  some  distance  removed  from  the 
exit  point  may  also  be  helpful  In  determining  the  merit  of  the  analysis. 

An  exarJ nation  of  the  exit  points  for  the  low  velocity  shots  (10-14)  shown 
xn  Fig.  2 reveals  that  they  are  approximately  4-6  inches  away  from  gages 
9-12  at  .e  center  of  the  plate.  The  maximum  strain  at  the  center  gages 


for  each  shot  is  2,5  x 10^  (10),  2,6  x 10^  (12),  3.3  x 10^  (13).  and 


3,2  X lO'^  (14)  roicrostraln*.  The  maximum  strain  from  the  SATANS  results  for 
shot  14  between  4 and  6 inches  from  the  exit  point  is  approximately  300 
microstrain  with  no  images  and  no  cavitation,  and  approximately  1000 


microstrain  when  cavitation  assuied.  Thus,  the  maximum  predicted  strain 
at  tl  a center  of  the  plate  is  underestifu^te  I, 

An  examination  of  all  of  the  measured  strain  traces  for  shots  10-14 

reveals  that  apparently  the  maximum  strain  for  all  five  shots  occurs  at 
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gage  4 and  is  approximately  -4.0  x 10  to  -4,5  x 10  microstrain**.  Also 
noted  is  the  fact  that  the  only  other  gage  that  registers  a large  negative 
(compression)  strain  for  these  shots  is  gage  These  two  gages  are  on  the 
wet  side  of  the  wall.  The  opposite  gages  on  the  dry  side  register  tension 
strains,  thus  indicating  large  bendin^  moments  at  these  two  gages.  This 
behavior  is  probably  due  to  the  fact  that  these  two  gages  are  closest  to 
the  clamped  edges  of  the  plate  where  the  bending  effects  would  be  pronounced. 
The  maximum  predicted  strains  near  the  edge  are  also  larger  than  those  at 
point  c, 

*The'  tracea  at  gages  9-12  were  loJt  for  shot  1_, 

**The  strain  traces  at  gage  4 look  essentially  the  same  as  that  of  Fig.  27 
^or  almost  every  shot,  including  the  high  and  medium  velocity  shots. 


As  a consequence  of  the  rather  poor  comparison  between  the  measured 
and  predicted  strains  when  cavitation  is  not  included,  and  the  extremely  long 
computation  time  required,  a full  scale  analysis  using  BR-IHR  has  not  been 
made.  The  results  from  preliminary  runs  using  BR-IHR  for  short  periods  of 
time  indicate  the  same  low  estimates  of  strain  that  SATANS  predicts. 


VII  SUMMARY  AND  CONCLUSIONS 


The  experimental  and  analytical  hydraulic  ram  and  structural  response 
programs  at  NWC  and  NFS  have  been  briefly  described.  Particular  attention  was 
given  to  the  hydraulic  ram  pressure  code  developed  by  Lundstrom  and  Fung  at 
NWC,  the  hydraulic  ram  tests  conducted  on  exit  walls  by  NWC,  the  computer 
study  of  the  exit  W2ill  response  by  NFS,  and  the  experimental  and  computer 
study  of  the  hydraulic  ram  shock  and  drag  phase  pressures  and  circular  entry 
wall  response  by  NFS.  The  goals  of  this  combined  NWC-NPS  effort  have  been 
the  development  of  an  understanding  of  the  hydraulic  ram-structural  response 
phenomenon  and  the  preparation  of  '‘.omputer  codes  that  can  be  used  to  accurately 
predict  the  fluid  pressures  and  the  subsequent  wall  response  to  the  penetrating 
projectiles. 

This  report  is  particularly  concerned  with  the  comparison  of  measured 
wall  strains  and  predicted  wall  strains  at  both  the  entry  wall  and  the  exit 
wall.  The  piston  theory  for  approximating  the  fluid-structure  interaction 
has  been  incorporated  into  the  two  computer  codes,  SATANS  and  BR-1,  allowing 
the  separation  of  the  fluid  and  structural  analyses.  For  the  entry  wall 
analysis,  the  Yurkovich  shock  pressure  theory  has  been  combined  with  the  NWC 
hydraulic  ram  pressure  code  to  give  a predicted  fluid  pressure  field  on  the 
the  entry  wall  as  a function  of  time.  This  predicted  pressure  field  was 
applied  to  the  circular  entry  wall  and  the  response  of  the  wall  was  computed 
using  SATANS.  A comparison  of  the  measured  and  computed  strain  histories 
twu  in.  from  the  center  of  the  plate  reveals  that  the  predicted  strains  have 
the  same  general  form  as  the  measured  strains,  but  the  magnitudes  of  the 
predicted  strains  are  smaller  than  the  measured  strains  shortly  after  the 
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response  begins.  When  cavitation  at  the  wall  Is  accounted  for  In  SATANS 
(no  net  tension  Interface  pressure) , the  magnitude  of  the  predicted  strains 
Increases  to  the  level  of  the  measured  strains.  However,  except  for  the 
circumferential  strain  on  the  dry  side,  the  similarity  In  the  fom  of  the 
strain  history  Is  lost.  In  addition,  the  predicted  radial  strain  with 
cavitation  shows  considerably  more  oscillation  than  the  measured  radial  strain 
Thus,  Its  questionable  whether  or  not  cavitation  occurs  at  the  entry  wall.  Th 
NWC  high-speed  motion  pictures  of  the  entry  wall  response  do  not  reveal  any 
cavitation  there. 

For  the  exit  wall  analysis,  the  predicted  Incident  fluid  pressure  field 
on  the  exit  wall  was  obtained  from  the  NWC  hydraulic  ram  pressure  code.  The 
Incident  pressure  was  assumed  to  drop  to  zero  when  the  bullet  exited  the 
tank.  The  other  walls  were  taken  either  as  free  surfaces  or  were  not  con- 
sidered In  two  separate  analyses.  The  response  of  the  20  Inch  square  exit 
wall  to  the  two  pressure  fields  was  computed  using  SATANS  by  Idealizing  the 
square  plate  as  a 20  Inch  diameter  circular  plate  centered  at  the  exit  point. 

A comparison  of  the  strains  In  the  exit  wall  measured  by  NWC  with  the 
strains  predicted  by  SATANS,  without  cavitation,  reveals  that  three  of  the 
four  predicted  strain  histories  have  essentially  the  same  form  as  those  of 
the  corresponding  measured  strains.  In  particular,  the  predicted  and  measured 
strain  behavior  appear  to  be  In  very  good  qualitative  agreement  at  the  time 
when  the  bullet  Is  estimated  to  exit  the  tank.  However,  the  magnitude  of  the 
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predicted  strains  Is  one  to  two  orders  of  magnitude  smaller  than  the  measured 
strains.  When  cavitation  Is  accounted  for  In  the  analysis,  the  magnitude  of 
the  predicted  strains  Increases  to  less  than  an  order  of  magnitude  smaller 


than  the  measured  strains,  but  the  similarity  in  the  form  of  the  strain 
histories  is  lost.  Thus,  it  appears  from  both  the  entry  wall  strain  com- 
parison  and  the  exit  wall  strain  comparison  that  the  cavitation  assumption 
of  no  net  tension  pressure  at  the  Interface  is  probably  not  correct.  A 
certain  level  of  tension  should  be  allowed. 

No  full  scale  analyses  of  the  exit  wall  were  conducted  using  the  modi- 
fied BR-1  code,  BR-lHP.,  due  to;  1)  a projected  minimum  execution  time  of 
about  3 hours  on  the  NFS  IBM  360/67,  2)  the  fact  t'  i.  rt  preliminary  runs 
using  BR-IHR  gave  the  same  level  of  small  strains,  as  and  3)  the  code 

in  its  present  form  cannot  handle  cavitation. 

As  a result  of  this  study,  the  major  conclusion  uun  be  made  that  the 
piston  theory  without  cavitation  yields  good  predictions  of  the  form  of  the 
strains,  but  the  magnitude  of  the  predicted  strains  is  often  an  order  of 
magnitude  or  more  too  small.  Addition  of  cavitation  of  the  fluid-structure 
Interface  leads  to  larger  predicted  strains,  but  the  similarity  in  the  form 
of  the  strains  is  lost.  Studies  are  currently  underway  to  modify  the  piston 
theory  to  allow  larger  strains  without  changing  the  form  of  the  strain 
history. 
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IX.  APPENDIX  - EVALUATION  OF  THE  BR-IHR  MODEL 


The  accuracy  of  the  BR-IHR  raodel  shown  In  Fig.  6 for  linear  response 
can  be  demonstrated  by  comparing  the  results  from  BR-IHR  with  the  series 
solution  to  the  problem  of  a simply  supported,  square,  fluid-backed  plate 
subjected  to  a uniform  step  pressure  load  P . The  piston  theory  equation 
governing  the  motion  of  the  plate  can  be  given  in  the  form 


where 


D it  + pew  ■ P 
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10.4  X 10®  psi 
0.0965  Ib/inJ 
1/3 
0.125 

1.0  lb/in.2 

20  in. 

386  in./sec^ 

7.64  X 10  #-8ec^/in. 


- Young’s  modulus 

- Specific  weight  of  the  plate 

- Poisson's  ratio 

- Thickness 

- Uniform  pressure 

- Plate  length  and  width 

- Acceleration  due  to  gravity 
Density  of  fuel 


*’•  « 53,000  in. /sec  - Acoustic  velocity  of  fuel 


The  solution  to  Eq.  (Al)  for  the  simply  supported  plate  can  be  given  in  the 
series  form 


w 


(A2) 


AWhen  executing  BR-IHR,  P was  taken  as  0.01  Ib/in.^  in  order  to  keep  the 
nonlinear  effects  small.  The  results  were  multiplied  by  100  in  the  compari- 
son with  the  series  solution. 
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where  the  top  sign  is  used  at  the  loaded  (wet)  surface  and  the  bottom  sign 
is  for  the  unloaded  (dry)  surface.  This  problem  is  selected  as  a test  case 
because  it  has  a relatively  simple  series  solution  and  it  roughly  approxi- 
mates the  actual  loading  condition  on  the  exit  wall. 

The  results  from  Eq.  A2  for  the  displacement  at  the  center  of  the  plate 
are  given  in  Fig.  A1  as  a function  of  time  for  pc  ■>  0 (no  fluid)  and  for 

pc  ■ 4.1  ^ . Note  the  significant  reduction  in  the  response  due  to  the 

in^ 

presence  of  the  fluid.  The  results  from  BR-lHR  for  the  displacement  at  the 

if  86  C 

center  of  the  plate  are  given  in  Figs.  A2  and  A3,  for  pc  " 0 and  4.1  — 

±n 

respectively,  for  50  time  steps.  The  time  increment  used  in  BR-IHR  was  the 
maximum  allowable  for  numerical  stability.  Also  shown  in  Figs.  A2  and  A3 
are  the  results  from  Eq.  A2.  Note  the  accurate  displacement  response  from 
BR-lIIR  for  both  values  of  pc. 

Profiles  of  the  plate  displacement  at  50  time  steps  with  pc  “ 0 are 
given  in  Figs.  A4a  and  A4b  for  y - 5 in.  and  x ■ 8 in.  respectively. 
Profiles  of  the  corresponding  outer  surface  strain  are  given  in  Figs. 

A5a  and  A5b  for  y - 6,5  in.  and  x » 10  in.  respectively.  Examination  of 
these  figures  reveals  that  the  displacements  are  accurately  predicted  by 
the  finite  element  model  but  the  strains  are  not.  In  particular,  the  strains 
in  the  center  portion  of  the  plate  are  Inaccurate.  However,  note  that  the 
strains  at  the  center  of  the  plate  are  relatively  small.  A plot  of  the  outer 
surface  strain  at  the  center  of  the  plate  as  a function  of  time  is  given 

in  Fig.  A6.  Note  the  relative  smoothness  of  the  finite  element  strain  com- 
pared to  the  strain  from  the  series  solution. 
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FIG.  A4.  PROFILE  OF  PLATE  DISPLACEMENT 


4 12 


20  28 


36 


44 


52  60 


(A5&)  x=IOin. 

FIG,  A5.  PROFILE  OF  PLATE  OUTERSURFACE  STRAIN 
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In  order  to  detennine  the  effect  of  the  model  upon  the  results,  a 144 
uiuie  model  was  also  considered.  A comparison  of  results  showed  that  the  144 
node  model  provided  somewhat  more  accurate  strains,  but  they  were  still 
inaccurate  in  the  center  portion  of  the  plate. 
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Wrlght-Patterson  AFB,  OH  45433 
Attn:  MAJ  John  W.  Mansur 
(AFFDL/FES) 

49.  J.  Rutland  (APATL/DLRV)  1 

Air  Force  Armament  Test  Lab 

Eglin  AFB,  FL  32542 

50.  Gerald  Bennett  (ASD/XRHD)  1 

Air  Force  Aero  Propulsion 

Laboratory 

Wrlght-Patterson  AFB,  OH  45433 

51.  Northrup  Corp.,  aircraft  Div,  1 

3901  Broadway 

Hawthorne,  CA  90250 
Attn:  Don  Wanlass,  F-18 
Code  2614/AH35 

52.  Columbus  Aircraft  Division  1 

Rockwell  International 

4300  East  Fifth  Avenue 
Columbus,  OH  43216 
Attn;  Robert  J.  Strohl 

53.  Andre  Holten  (AFFDL/FES)  1 

Air  Force  Flight  Dynamics 

Laboratory  (AFSC) 

Wrlght-Patterson  AFB,  CH  45433 

54.  Dr.  M.  J,  Jacobson  1 

Orgn  3761/62 

Northrop  Corp. 

Aircraft  Division 
3901  Broadway 
Hawthorne,  CA  90250 
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